Local Wettability Reversal during Steady-State Two-Phase Flow in Porous Media 
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We study the effect of local wettability reversal on remobilizing immobile fluid clusters in steady- 
state two-phase flow in porous media. We consider a two dimensional network model for porous 
medium and introduce a wettability alteration mechanism. A qualitative change in the steady 
state flow patterns, destabilizing the percolating and trapped clusters, is observed as the system 
wettability is varied. When capillary forces are strong a finite wettability alteration is necessary 
to move the system from single-phase to two-phase flow regime. For the case of both phases being 
mobile we find a linear relationship between fractional fiow and wettability alteration. 

PACS numbers: 47.56. -(-r,47.61.Jd 



In a time of dwindling oil reserves, the fact that some 
20 to 60 percent of the oil remains unrecovered at the end 
of oil production constitutes a challenge of increasing im- 
portance Jj]. The reason for this loss is the formation of 
oil clusters embedded in water and held in place by cap- 
illary forces, which in turn are controlled by the wetting 
properties of the reservoir fluids with respect to the ma- 
trix rock. Therefore, the wetting properties of crude oil, 
brine and rock systems is a central research topic in the 
petroleum industry within the field of Enhanced Oil Re- 
covery (EOR). EOR covers recovery methods other than 
pressure depletion (typical primary recovery) and water 
or gas injection (typical secondary recovery). 

Formation wettability was reviewed recently in 
Reservoir rocks do not have uniform wetting properties. 
When a reservoir fluctuates spatially between strongly 
oil-wet and strongly water-wet it is called fractionally wet 
Sandstone is strongly water- wet before oil migrates 
from a source rock into the reservoir. After oil enters a 
pore it forms water films 0, and a wettability alter- 
ation of the rock is believed to take place which can be 
irreversible due to direct adsorption of asphaltenes d, 0] 
or can be reversible due to other surface active molecules 
present in the crude oil, such as caboxyl acids. 

Both laboratory experiments and field tests have 
shown that deviation from strongly oil-wet to water-wet 
or neutral-wet conditions significantly increase oil recov- 
ery efficiency 0. Wettability of reservoir rocks can be 
altered by changing the brine composition, e.g., lower- 
ing safinity adding water-soluble surfactants j[lo| or 
even by adding oil-soluble organic acids or base [8|. An 
increase in temperature also increases water wetness of 



the reservoir [llj. Correlations have been shown with 



wetting behavior to the electrostatic forces between the 
mineral and oil surfaces [l2| . Increasing the salinity level 
of the pore water makes the wetting transition shift in 
pH. The oil industry therefore paying a great attention 
towards low-salinity water flooding. 

So far, there is no consensus on the dominating micro- 
scopic mechanisms behind the transport of the agent that 



mediates the change in the pore wettability. The result 
of it, however, will be that there are correlations between 
the positions of the altered pores. In light of the mas- 
sive uncertainty connected to the transport mechanisms, 
we simplify the problem assuming no correlation between 
positions of the altered pores. The precise question we 
pose in this paper is the following: Will the reversal of 
the wetting properties of a percentage of the pores lead to 
a significant change in the fractional flow patterns? We 
will see here that even without wettability switch cor- 
relations there can be a permanent change in fractional 
flow when the switched fraction of pores is large enough. 
How spatial correlations between the altered pores fur- 
ther affect the flow is a matter for future investigations. 

In order to shed light on how local wettability alter- 
ations may remobilizc stuck fluid clusters and to look for 
the consequences of this on the flow properties, we study 
here the effect of such alterations on a two-dimensional 
pore scale transport network model [l3l - [l5| . The mech- 
anism for local wettability alteration is introduced in 
the model by changing the wetting contact angles of the 
menisci at pore necks. To study steady-state flow reflect- 
ing the situation deep inside the reservoir, we implement 
biperiodic boundary conditions. The network is repre- 
sented by a square lattice of cylindrical tubes, tilted 45° 
with respect to the imposed pressure gradient as well as 
the overall flow direction. Disorder is incorporated in the 
system by assigning the radius (r) of each tube randomly 
from a flat distribution in the range [O.IZ, 0.4^], where I is 
the length of the tube. The network is filled with two im- 
miscible fluids, named as water and oil, which flow inside 
the tubes (links). With respect to a pore, one fluid is less 
wetting while the other is more — hence there is a con- 
tact angle 9 associated with each meniscus. The tubes 
are cylindrical with respect to permeability, but they are 
considered hour-glass shaped with respect to the capil- 
lary pressure Pc- At position x, Pc is obtained from a 
modified form of the Young-Laplace law 
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where 7 is the interfacial tension between the fluids. We 
define the wettabihty of the system by w = -/Voii/iVtotai, 
where A'oii is the number of oil- wet hnks among the total 
iVtotai links in the system. Thus lo — and lo — 1 rep- 
resent pure water-wet and pure oil-wet systems respec- 
tively. A fractional value of w represents a fractionally 
wet system. When the wettability of a pore is altered, 
the direction of capillary forces at the menisci inside that 
pore changes. We therefore change the wettability of a 
pore in the model by altering the sign of Pc- In this work 
we consider a mixture of perfectly oil-wet and perfectly 
water- wet conditions, i.e., 9 is either 0° or 180°. 

The flow is driven by setting up an external global 
pressure drop. The local flow rate g in a tube with a 
pressure difference Ap between the two ends of that tube 
follows the Washburn equation of capillary flow 1^ 12 1; 
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where k is the permeability. For cylindrical tubes, k — 
which is known from Hagen-Poiseuille flow. Here 
a is the cross-sectional area of the tube and ^ is the 
volume average of the viscosities of the phases present 
inside the tube. The sum over pc runs over all menisci 
within the tube. The net flow at a node must be zero 
in a network. Eq. [2] then gives a large number of linear 
equations connecting the pressure variables at the nodes. 
The pressure field is obtained by solving the linear equa- 
tions and the flow fleld then follows from eq. [2J The 
system is forward time integrated by one step, using the 
explicit Euler scheme. Inside a tube all menisci move 
with a front speed determined from the local flow rate 
q, but when a meniscus reaches the end of a tube new 
menisci are formed in neighboring tubes maintaining the 
volume conservation. A maximum number of six menisci 
is considered here in a given tube. Further details of how 
the menisci are moved can be found in [l4| . though they 
are not necessary to understand the results of wettability 
alteration here. 

The flow is controlled by the ratio between capillary 
and viscous forces at the pore level and quantifled by 
capillary number Ca — /i(5tot/(7S) where S is the cross- 
sectional area of the network and Qtot is the total flow 
rate. Two different capillary numbers Ca = 1.0 x 10~^ 
and Ca = 1.0 x lO^'^ are considered here. The ratio be- 
tween the viscosities of the two fluids, M = /Iq/Mw is an- 
other dimensionless parameter to control the flow, where 
Ho and /iw are the viscosities of oil and water respectively. 
Here we restrain our discussion to viscosity matched flu- 
ids which gives M = 1. Simulation is done considering a 
network of 40 x 40 links. This is sufficiently large to be 
in the asymptotic limit for the range of parameters 14 1 . 
An average over 10 different samples has been taken for 
each simulation. 

The simulation is performed under constant flow rate 
Qtot which sets the capillary number. Initially, the sys- 





FIG. 1. Distribution of normalized local flow rates (q) (a) 
before and (b) after the wettability alteration at Ca — 1.0 x 
10"^ and S = 0.6. The long arrow indicates the overall flow 
direction. The highest to lowest flow rates are indicated by 
darkest to lightest gray as shown by the scale. The steady 
state flow before wettability alteration is dominated inside 
only a few percolating paths as seen in (a). In (b), wettability 
of 50% tubes are altered and the percolating path breaks, 
causing the flow to be distributed over the whole network. 



tem is pure water- wet (w = 0) and filled with given oil 
and water saturations which remain constant throughout 
the simulation. Thus the saturation, capillary number 
and the wettability are the independent variables here. 
As the system evolves with time, it either starts a viscous 
fingering for high values of Ca or invasion percolation pro- 
cess when capillary forces dominate. Due to biperiodic 
boundary conditions, both drainage and imbibition take 
place simultaneously, leaving wetting and non-wetting 
fiuid clusters in the system. On long timescales the sys- 
tem evolves to a steady state characterized by the sys- 
tem's macroscopic properties, such as global pressure, 
that becomes constant on the average. Mixing of two 
phases in steady state depends on the external control 
parameters, both Ca and the saturation of fiuids. In 
Fig. [Ija), the distribution of normalized local flow rates 
at steady state is shown for the lower capillary number 
Ca = 1.0 x 10""^. If the difference in the saturation of the 
two fluids is large, one of the phases percolates the sys- 
tem with trapped immobile clusters of the other phase as 
can be seen in Fig. [Ua). Clearly, the flow is dominated 
through a few anisotropic percolating paths [l^ : leaving 
trapped immobile clusters in the system. In this case, we 
may apply Darcy's law for the single-phase pressure 
along the system length L by. 
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which will be significantly larger than that of single-phase 
flow due to the smaller effective cross section E'. How- 
ever, in the other situation when neither phase percolates 
the system and both the phases are mobile, all spanning 
flow paths contain interfaces associated with capillary 
forces. This situation requires even larger driving pres- 
sure due to the local dynamics of interfaces and capillary 
barriers. Therefore, as soon as low resistance percolating 
path appears in the system the global pressure drop de- 
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FIG. 2. Variation of oil fractional flow with time. Three 
different regions in each curve imply the time evolution in 
a water-wet system, then in a fractionally wet system after a 
wettability alteration cj and finally restoring the initial water- 
wet system. For each saturation 10 different curves from top 
to bottom correspond to different values of u} from 0.1 to 1.0 in 
intervals of 0.1. As u> increases a gradual change in fractional 
flow is observed. 



creases and the percolating and trapped clusters remain 
static. This is not favorable for oil-recovery purposes as 
it leaves immobile fluid in the reservoir. 

After the system evolves to a steady state we change 
the vifcttability w by altering the direction of capillary 
forces in required number of links, chosen randomly. This 
initially perturbs the pressure field, deviates the system 
from existing steady state and subsequently leads it to 
settle into a new steady state. This again results to break 
any existing percolating and trapped fluid clusters which 
start flowing within a few time steps as illustrated in 
Fig.[TJb). If the wettability alteration is sufficiently large 
the change to a non-percolating system is permanent, 
otherwise new percolating cluster eventually appears. 

We measure the fractional flow of oil, defined as the ra- 
tio between the flow of oil (Qoii) with the total flow (Qtot) 
as F^{S) = Qow/Qtot in a system with wettability w and 
oil saturation S. Independent simulations have been per- 
formed for different values of uj and S. In Fig. [51 F'^{S) 
is plotted as a function of time. The three regions of the 
curves correspond to the initial pure water- wet {uj — 0) 
system, flipping wettability by a nonzero value of uj and 
then again flipping them back to the initial lu — value. 
There is a few things to notice. First, there is a grad- 
ual decrease in the steady state oil fractional flow with 
increasing w, i.e., more oil-wetness of the system. This 
is consistent with experimental results as it is observed 
that changing from oil-wet to mixed and water-wet sys- 
tems, as in the case for low-salinity flooding, increases oil 
production [l^] ■ Secondly, the alteration process is re- 
versible here, as in the third regime of the curves F'^{S) 
returns to the initial steady state value F^{S). Third, 
when the capillary forces are strong and the system is in 
a single-phase flow regime, a finite amount of wettabil- 
ity alteration is needed to move it towards a two-phase 
regime. This can be seen in Fig. [Hb), that for low cap- 
illary number Ca = 1.0 x 10~^ and saturation S — 0.7, 
apart from an initial transient response, the steady state 



FIG. 3. Fractional flow as a function of saturation in water- 
wet (O) and oil-wet (□) systems. 1 — ^"(1 — S), calculated 
from F'^IS), is drawn by the dashed line which matches nicely 
with the F'^{S) curve. 



fractional flow does not change for lu < 0.2. 

The variation of oil fractional flow with oil saturation 
S" in a pure water-wet and a pure oil-wet systems are 
plotted in Fig.[3l In both the systems, the fractional flow 
does not follow F = S value represented by the solid 
straight line. This is due to the presence of interfaces 
with capillary pressure. The curves cross the F — S line 
from below at some point less than 50% saturation in 
pure water-wet and higher than 50% saturation in pure 
oil-wet systems. Fractional flow always higher in water- 
wet system {F^{S)) than in the oil-wet system {F^{Sj). 
The difference AFi(S') = F°{S)-F^{S) depends on the 
saturation, it is zero at the two end points S — and 
1. It is maximum in between, where both the phases 
are mobile. The maximum value of AF increases with 
decreasing Ca, it is « 0.22 for Ca = 1.0 x lO'^ and w 0.40 
for Ca — 1.0 X 10~'^. This due to the shape of the curves 
which strongly depend on Ca 2^, 21 1. 

As we consider only the case of viscosity ratio M — 
1, a symmetric relationship between the fractional flow 
in pure water-wet and pure oil-wet systems is expected. 
From straight-forward symmetry argument we write the 
relation between the oil fractional flow in the two different 
systems as. 



F\S) = 1~F°{1~S) 



(4) 



In Fig. H 1 - F°(l - S") is plotted by the dotted hues and 
it matches nicely with F^ (S) curve. 

The wettability of real oil reservoirs generally have 
fractional wetting conditions. We therefore measure 




FIG. 4. Change in oil fractional flow A_F"(S) with saturation 
and system wettability. For the lower capillary number, a 
finite UJ is required to have a non-zero value of A_F'"(5') at 
high saturations. 
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FIG. 5. Change in the oil fractional flow A_F'"(5') as a function 
of uj at intermediate range of saturation, from S — 0.3 to 0.7. 
It can be seen that AF^{S) varies linearly with to. 



F'^{S) in different systems with fractional values of uj 
from 0.1 to 1.0 in intervals of 0.1. The difference of 
F'^{S) from that in pure water- wet system, AF'^{S) = 
F"{S) - F'^{S), is then calculated. In Fig. H AF'^{S) 
is plotted as a function of saturation and wettability for 
the two different capillary numbers. It can be seen that 
the surface is more flat for the higher capillary number. 
Moreover, for the lower Ca there is a finite region with 
AF'^{S) = to be noticed. It is already observed in Fig. 
[Sfb) that for high saturation a finite wettability alter- 
ation ujtr is required to permanently change the system 
from a single-phase dominated regime. For w < Lu-tr we 
have AF'^iS) = 0. 

In order to find the functional dependence of AF'^{S) 
on UJ in the regime where both the phases are mobile, 
AF'^{S) is plotted against w in Fig. [5] for saturations 
from S = 0.3 to 0.7. The variation of Ai^"(5) with lu is 
found to be linear here with some fluctuations. For w — 0, 
AF°{S) =0. For = 1, AF^{S) directly follows from 
Fig. [las the difference between the F°{S) and F'^{S) 
curves. Assuming linearity in between, one can write 

AF'^ (S) = oj AF\S) = u}[F°{S) - F\S)], (5) 

which we can combine with eq. |4] to obtain 

F'^iS) = F°{S)-io[F"{S)+F"{l -S)-l]. (6) 

This gives the functional dependence of fractional flow 
on both saturation and system wettability. It connects 
fractional flow in a fractionally wet system with that of a 
purely water-wet system by a linear relation with uj and 
S. As fractional flow curves may be more easily obtained 
in samples of pure wetting properties, this equation may 
be useful to have an estimate of fractional flow when 
wettability is altered by, e.g., low-salinity water injection. 

As AF^{S) is never negative it is clear that the frac- 
tional flow of oil will never decrease with the increase of 
water-wetness of the system. For a large viscosity con- 
trast this result could vary, as our symmetry argument 
requires viscosity matched fluids. It is therefore requires 
further extensive numerical simulations to identify the 
results in that regime. 



As a conclusion, we study the effect as well as the de- 
tailed pore level mechanism that causes the change in the 
steady state two-phase flow properties due to wettabil- 
ity alteration in porous media. A change in the system 
wettability causes a perturbation in the system's flow 
pattern to destabilize any percolating and trapped im- 
mobile clusters appeared in the steady state. In order 
to prevent forming similar structures again, a sufficiently 
strong wettability alteration is required depending upon 
the capillary number and saturation. When both the 
phases are mobile, fractional flow of oil is found to in- 
crease linearly with increasing water- wetness of the sys- 
tem. The results of our simulation are general in nature 
and consistent with the experimental results reported in 
literature qualitatively. They show that the remobiliza- 
tion of fluid clusters is due to purely fluid mechanical 
changes in the system induced by changes in the wetting 
contact angle — which in turn comes from the wettability 
alterations of the porous medium. 
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